Ultracold atoms in optical lattices with random on-site interactions 
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We consider the physics of lattice bosons affected by disordered on-site interparticle interactions. 
Characteristic qualitative changes in the zero temperature phase diagram are observed when com- 
pared to the case of randomness in the chemical potential. The Mott-insulating regions shrink and 
eventually vanish for any finite disorder strength beyond a sufficiently large filling factor. Further- 
more, at low values of the chemical potential both the superfluid and Mott insulator are stable 
towards formation of a Bose glass leading to a possibly non-trivial tricritical point. We discuss fea- 
sible experimental realizations of our scenario in the context of ultracold atoms on optical lattices. 

PACS numbers: 03.75.Ss,05.30.Fk 



X 



Ultracold atomic gases have attracted large interest 
during the last years, in particular due to the experimen- 
tal achievement ofquantum degeneracy both for bosonic 
P| and fermionic Q gases. Very recently, the extraordi- 
nary experimental developments in the manipulation of 
ultracold atoms have opened the path towards a new fas- 
cinating research area, namely the analysis of strongly- 
correlated atomic gases. In this sense, remarkable ex- 
perimental results have been already reported, including 
experiments on the so-called BEC-BCS crossover in de- 
generate Fermi gases and experiments on ultracold 
gases in optical lattices, such as e.g. the observation of 
the Mott insulator-to-superfluid transition 0, [j| and the 
experimental realization of a Tonks-Girardeau gas [|J . 

Up to now the experiments performed in optical lat- 
tices have considered a defect-free square lattice poten- 
tial. However, recent proposals have discussed new sce- 
narios which depart from this standard. In particu- 
lar, various lattice geometries are achievable by stan- 
dard laser techniques, and should lead to fascinating 
physics 0,0- On the other hand, disorder may be in- 
duced in the laser potential in different ways. Local- 
ized impurities can be created, leading to Kondo-likc 
physics 0| . Randomness can be produced by means of 
additional incommensurable lattices, or by lase r sp eck- 
les [ljj, an( i ma y l ea( i to Anderson localization [ll[ and 
Bose glass phases [l^. First experimental steps towards 
random potentials have already been achieved Ad- 
ditionally, randomness may be induced in the hopping 
rates in the optical lattice, leadin g to Fermi glasses, spin 
glasses and quantum percolation [14|. 

One of the most fascinating possibilities for the con- 
trol of atomic gases is provided by the manipulation of 
the interatomic interactions by means of Feshbach reso- 
nances jl^- Using a suitably adjusted magnetic field a 
resonance is induced between a molecular state and the 
unbound states of the incoming atoms, leading to im- 
portant modifications of the scattering properties. The 
s-wave scattering length undergoes a very large change 
in a relatively narrow window of values of the applied 
magnetic field, becoming ±oo at the resonance. This ex- 



treme sensitivity of the scattering properties at the verge 
of a Feshbach resonance is crucial for our discussion. 

In this Letter, we consider the physics of lattice bosons 
subject to a novel kind of disorder, namely bounded dis- 
order in the strength of the interatomic interactions |ln1 | . 
In particular, we show how this type of disordered system 
can be realized in the context of ultracold atoms near a 
Feshbach resonance. We discuss the consequences that 
the disorder in the interaction strength has on the zero- 
temperature phase diagram of the system by means of 
a strong-coupling expansion (SCE) and quantum Monte 
Carlo (QMC) simulations, and contrast our findings to 
the case of randomness in the chemical potential 12]. 

The above scenario can be realized in a gas of ultracold 
bosons confined to an optical lattice using state-of-the- 
art experimental techniques on atom chips ^j- I n the 
following, we consider a one-dimensional configuration, 
where the atoms are assumed to be strongly confined 
in the other two dimensions (by e.g. additional optical 
confinement). However, qualitatively similar results as 
those discussed in this Letter also hold in higher dimen- 
sions. We assume that the gas is brought at the verge of 
a Feshbach resonance by an off-set magnetic field, where, 
as discussed previously, slight modifications of the mag- 
netic field lead to large variations of the scattering prop- 
erties. We furthermore consider the bosonic gas to be 
close to a magnetic wire, inducing a spatially random 
magnetic field |l7j. which can be considered sufficiently 
weak, such that the variations in the potential energy can 
be considered negligible when compared to other energy 
scales. However, since the off-set field sets the system 
at the verge of a Feshbach resonance, the slight random 
variations of the magnetic field lead to a spatially ran- 
dom variation of the local interatomic interactions (see 
Fig. ^) . We assume that the variations of the scattering 
properties are bounded, i.e. the spatial variations of the 
magnetic field do not induce a crossing of the resonance. 

In order to avoid significant van der Waals losses, the 
atoms should be placed at a distance > 0.5/im from the 
surface of the atom chip, which restricts the character- 
istic wavelength of the variations of the magnetic field 
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FIG. 1: Typical spatial variation of the on-site interactions. 
In the figure we have considered a 87 Rb cloud 0.5/im over the 
surface of an atom chip, at a 1/im-wide magnetic wire with a 
current leading to a magnetic field variation of up to ±2.5mG, 
an off-set field 5mG away from the 685. 8G resonance, a lattice 
spacing of lfi, and a lattice depth of 5 recoil energies. Ui is 
obtained by integrating J dxU(x)\ipi(x)\ 4 / J dx\ipi(x)\ 4 , with 
ij}i(x) the i-th on-site wavefunction which is approximated by 
a Gaussian. 

at the wire to Ax > l/_tm. The inter-site separation in 
the optical lattice is approximately A = 400nm, but can 
be increased to 1/zm by setting an angle between the 
counter-propagating lasers in such a way that Ax < A. 
The variation of the magnetic field at the wire can be 
adapted to typical Fcshbach resonance widths, which are 
of the order of 1 — 10 mG. Fig. shows an example of 
the variation of the on-site interaction for a typical ex- 
perimental setup. 

Under these conditions, the system is described by a 
Bose-Hubbard Hamiltonian with random interatomic in- 
teractions: 

<ij> i i 

(1) 

where bi (b\) is the annihilation (creation) operator for 
bosons in the i-th lattice site, = b\bi, J is the hopping 
constant, jii = [i is the local chemical potential, and 
Ui denotes the local value of the coupling constant for 
the interatomic interactions. We assume U uniformly 
distributed inside the range U(l — e) < [7j < U(l + e) 

m 

In order to assess the effects of randomness in the 
on-site interactions and to contrast it to the case of a 
random chemical potential, we shortly review the latter 
case, where \i — A < m < \i + A, at a fixed U% = U 
(e = 0) 0]: The phase diagram in the strongly inter- 
acting regime can be obtained from considering first the 
trivial limit, J = 0: One finds that if fi falls within an 
interval (n — 1)U + A < fj, < nU — A, with an integer n, 
the ground-state of the system is an incompressible Mott 
insulator with n bosons on each lattice site. For values of 
[i outside these intervals, i.e. for nU — A < [i < nU + A, 
the system enters an insulating, but compressible phase 
commonly referred to as Bose glass, in which the occupa- 



tion per site varies between n and n + 1. This physics is 
maintained up to moderate values of J, leading to a phase 
diagram with three distinguished phases: superfluid (SF) 
at large J, Bose glass (BG), and the Mott insulator (MI) 
regions. Recent quantum Monte Carlo calculations have 
shown that the transition from SF to MI always occurs 
through an intermediate BG phase 0], as suggested in 
Ref. O. 

A characteristic difference to the case of random U at 
a fixed value of /ij = fi (A = 0), is already obtained 
in the limit J — 0: Now compressible phases occur for 
n(l — e) < iijU < n(l + e), in which the occupation per 
site varies between n and n + 1, whereas outside these 
regions the system enters into a MI. In contrast to the 
case of randomness in the chemical potential, the extent 
of compressible regions between the MI thus increases 
with the occupation number n. The disorder in U thus 
leads to an increased destabilization of MI states with 
increasing n. Eventually, MI regions with an occupation 
of n > (1 + e)/2e particles per site disappear. There- 
fore, given a disorder strength e, only a finite number of 
MI regions remain stable. This selective destruction of 
MI states cannot be accessed in the case of a random 
chemical potential, where all the MI lobes vanish once 
A > U/2. More generally, if both U and fii have a 
bounded random distribution, the n-th lobe disappears 
for (2n- l)e + 2A/U= 1. 

In order to analyze the phase diagram for finite hop- 
ping J > 0, we first estimate the extent of the MI regions 
using the SCE 20] . This method allows for a quantita- 
tive calculation of the boundaries between compressible 
and incompressible phases for sufficiently low J, particu- 
larly in low dimensions, where the coordination number 
is small. 

From the SCE one obtains the energy gap for adding 
or removing a particle from the MI, and performs an 
expansion in the small parameters J/U, A/U and e. For 
a given value of J, these gaps determine the boundaries 
of the MI lobes, which are obtained in first order as 

£ = -2(n+l)^-£+n(l- C ) (2) 

§ - 2n^ + ^ + (n-l)(l + e) (3) 

for the upper (u) and lower (I) boundaries of the MI with 
integer filling n. In particular we find, that the lower 
boundary of the first MI lobe (n = 1) is not affected by 
the disorder in U, in contrast to the random- /x case |2(1. 
This reflects the observation, that no compressible phase 
emerges below the n — 1 MI at J = for the random- U 
case, in contrast to the case of random fj,. 

From equating the two boundaries one obtains an 
(under-) estimate for the largest extent of the MI phase 
with filling n: 

J_c_ = 1 - 2A/U - (2n - l)e 

U 2(2n + l) ' ' 
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Calculating higher-order terms in the strong-coupling 
expansion improves the quantitative description of the 
MI lobes, but does not affect our qualitative conclusions. 
In particular, all the higher order terms in the random 
interaction are at least of order (J/U) 2 e. Considering 
that the critical hopping J c /U <C 1, even relatively large 
disorder strengths may be approximated using the first 
order shift of the phase boundaries. For the figures pre- 
sented below, we used expansions up to third order and 
also considered finite size effects |20j. The latter, as 
shown below, are rather significant, and converge only 
very slowly to the thermodynamic limit when the num- 
ber of lattice sites considered increases. 

While the SCE allows to estimate the extent of the 
Mott-insulating phases quantitatively, it does not provide 
information about the complete phase diagram, in par- 
ticular concerning the nature of the compressible phases 
and the presence of SF and BG phases. In order to obtain 
a more complete phase diagram for bosons with random 
on-site interactions we have performed QMC simulations 
using the stochastic series expansion method |2j| with di- 
rected loop updates [22I . We performed simulations 
for periodic chains with L = 200 sites averaged over typi- 
cally 200 disorder realizations for each data point. In the 
simulations, the temperature T was chosen low enou gh to 
obtain ground state properties of the finite systems [2J| . 
In order to reliably identify the various phases of the sys- 
tem, we measured the global compressibility k = dN / dp 
from the fluctuations of the particle number N, and the 
winding number fluctuations (W 2 ), from which the su- 
pcrfluid density is obtained as ps = LT/(2J) (W 2 ) 
Fig. [3 shows the disorder averaged values of both ob- 
servables as functions of J/U along a cut of constant 
p/U = 0.51 for e = 0.25, which at J = belongs to the 
MI region with n = 1. At low values of J, an extended 
MI region with n = 1 is clearly identified by vanishing 
values of both k and p$- For large value of J/U, both k 
and ps take on finite values, thus identifying the large- 
J region as a SF phase. In the intermediate region, for 
0.078 < J/U < 0.133, the system shows a finite com- 
pressibility, but no superfluid response. This regime is 
thus characterized as a disorder-induced BG phase, sep- 
arating the MI and SF region at this value of \x/U = 0.51. 

Performing similar scans along different lines in pa- 
rameter space, we obtain QMC estimates of the various 
phase boundaries, leading to the phase diagram of Fig. 
Here, we combined the QMC results for e = 0.25 with 
third-order SCE results for the extent of the MI phases 
in the thermodynamic limit (TDL). Within the chosen 
resolution, the finite size corrections of the phase bound- 
aries obtained by QMC show only a weak dependence on 
the system size. An accurate determination of the MI 
boundaries from QMC in the disordered case [If| would 
require significantly larger s yste m sizes due to rare re- 
gions of delocalized bosons 20] . In fact, performing a 
third-order SCE for a finite chain of L = 200 sites com- 
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FIG. 2: Zero temperature compressibility k and superfluid 
density ps of bosons on a one-dimensional optical lattice with 
random interaction strength of e = 0.25 as a function of J/U 
along a cut a p/U = 0.51, obtained from QMC simulations. 
The error bars quantify the standard deviation of the sample- 
to-sample fluctuations, whereas the dashed lines indicate the 
estimated extent of the Bose glass along this cut. 
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FIG. 3: Zero-temperature phase diagram of bosons on a one- 
dimensional optical lattice with random interaction strength 
of e = 0.25, obtained from QMC simulations of 200 sites and 
a third-order SCE for the thermodynamic limit (TDL). The 
extent of the Mott lobes in the pure case (e = 0) from SCE 
is indicated by the dashed line, whereas the dot-dashed line 
shows the SCE results for a finite system of L = 200 sites. 

pares well to the estimated extend of the MI regions from 
QMC (Fig. in the strong-coupling regime, J/U < 0.1, 
and indicates the rather strong finite size effects in the 
disordered case. 

Our numerical calculations confirm the qualitative pic- 
ture considered above: The MI lobes clearly shrink, not 
only in the /x-direction but also in J, and eventually dis- 
appear for sufficiently large values of p. In the case con- 
sidered in Fig. |3 where e = 0.25, only the MI lobes 
with n = 1 and 2 remain stable. Furthermore, we ob- 
serve that the lower boundary of the n = 1 MI does not 
vary significantly with the disorder strength up to the 
critical hopping J c . The reduced relevance of disorder 
in the interaction strength in the low-density region is 
also reflected by the complete absence of a BG phase for 
p < 0. Therefore, a dilute compressible lattice boson gas 
remains superfluid even in the presence of disordered in- 
teractions. This is in clear contrast to the case of random 
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\i, where the BG extends down to arbitrary low but fi- 
nite densities. The absence of a BG for [i < indicates a 
tricritical point for the SF, MI, and BG phases along the 
lower boundary of the n = 1 MI. However, our calcula- 
tions cannot reliably distinguish whether the BG phase 
disappears exactly at a tricritical point at fj, = J = 0, or 
if this occurs for finite values of /x, J > 0. This interesting 
question thus remains open for future studies. 

To summarize, in this Letter we analyzed the novel 
physics occurring when the on-site interaction in a Bosc- 
Hubbard Hamiltonian acquires a bounded random char- 
acter. We discussed how this exciting possibility can be 
realized in ultracold atomic gases in optical lattices at the 
verge of a Feshbach resonance, when they are brought 
in the vicinity of weak spatially random magnetic fields 
generated by magnetic wires. As for the case of ran- 
dom chemical potentials, a phase diagram with super- 
fluid, Bose-glass and Mott-insulator phases is predicted. 
However, important differences can be found when com- 
paring the random U and random /x cases: For random 
interactions, the Mott-insulator phases become progres- 
sively narrower (both in the \x and J direction) for larger 
occupation numbers, and eventually disappear beyond 
a given filling. Furthermore, the Bose-glass phase disap- 
pears for low chemical potentials and hopping, and hence 
a tricritical point occurs, although our numerical calcu- 
lations cannot discern whether this occurs at fj, = J = 
or very small but finite values. 

Our calculations were performed in absence of an over- 
all (usually harmonic) confining potential, whereas the 
necessary commensurability for Mott-insulator phases 
can be achieved in practice only in the presence of in- 
homogeneous potentials. However, if the confining po- 
tential is sufficiently shallow, a local chemical potential 
may be considered, and local phase occur, as tho roug hly 
demonstrated by means of QMC calculations [2jJ 12£| for 
the pure case. Of course, the finite size of these differ- 
ent regions could significantly affect the boundaries be- 
tween the Bose glass and supernuid, once the localization 
length becomes smaller than the spatial extent of the var- 
ious regions. Typical experiments on the Mott-insulator 
phase of atoms in optical lattices 0, 0| rely on the obser- 
vation of the broadening of interference fringes in time- 
of-flight pictures (provided by the insulating character 
of the phase), and the monitoring of the opening of the 
Mott-gap in the excitation spectrum by means of tilt- 
ing experiments. In these experiments, although other 
spatial phases are expected, the Mott-insulator region is 
the largest one and dominates the experimental detec- 
tion. Similar experiments could be performed to detect 
the shrinking of the Mott-insulator phases as a function 
of the disorder in U predicted in this paper. 
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